Interplay of coupling and superradiant emission in the optical response of a double 

quantum dot 
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We study theoretically the optical response of a double quantum dot structure to an ultrafast 
optical excitation. We show that the interplay of a specific type of coupling between the dots and 
their collective interaction with the radiative environment leads to very characteristic features in the 
time-resolved luminescence as well as in the absorption spectrum of the system. For a sufficiently 
strong coupling, these effects survive even if the transition energy mismatch between the two dots 
exceeds by far the emission line width. 
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I. INTRODUCTION 

Systems composed of two quantum dots (QDs) have 
attracted much attention in recent years. Many theoreti- 
cal and experimental results have demonstrated that the 
physical properties of such double quantum dots (DQDs) 
are much richer than those of individual ones. On the 
one hand, this may pave the way to new applications, in- 
cluding long-time storage of quantum information^, con- 
ditional optical control 2 that may lead to an implemen- 
tation of a two-qubit quantum gate 3 -, generation of en- 
tangled photons 4 or coherent optical spin control and 
entangling-^ 7 -. On the other hand, in order to take ad- 
vantage of these extended possibilities, the properties of 
DQDs have to be understood and controlled. 

There are two major factors that determine the physi- 
cal (in particular, optical) properties of DQDs: the cou- 
pling between the dots and their interaction with the en- 
vironment. In both these areas new features appear, as 
compared to the physics of individual dots. Both prob- 
lems have been in the focus of extensive experimental 
and theoretical work but many important questions re- 
main open. 

Coupling between the dots is essential for quantum 
conditional control, entanglement formation and imple- 
mentation of two-qubit gates. It is therefore under- 
standable that considerable experimental effort has been 
devoted to demonstrating its existence in double dot 
systems 8 -! 9 - ' 10 ' 11 ' 12 , while theoretical studies were aimed 
at characterizing its signatures in the optical response of 
DQDsI 3 -. 

Coupling of DQDs to their environment is affected 
by collective effects, which may lead to superradiance 
phenomena 1 ^. Theoretical studies on the dephasing (in 
particular, decay of entanglement) in double dot systems 
have shown that coherence properties strongly depend on 
whether the dots are coupled to a common reservoir or 
to separate reservoir a 15 ' 16 ' 17 ' 18 . The collective nature of 
the coupling to the environment allows one to construct 
arrays (collective quantum bits) which are more resistant 
to decoherence than individual system o 19 ' 20 . In fact, in 
any such array of two-level systems collectively coupled 
to their common reservoir, the dephasing of some states 



is slowed down, while other states decohere faster. In the 
case of radiative decay, an essential role is played by the 
transition energy mismatch between the systems forming 
the array^i. These subradiance and superradiance effects 
influence the optical response of DQDs and can affect the 
coherence properties of DQD-based quantum devices. 

Since the paper by Dicke^, coherent effects in the ra- 
diative decay of two or more atoms have been studied in 
numerous works. The emission from identical— i 24 i 25 i 26 i 27 
and non-identica l 28 ' 29 two-level systems has been studied 
and methods suitable for the description of arrays of var- 
ious shapes have been develope d 30 ' 31 . Systems formed 
by QDs share many features with those made of natu- 
ral atoms. In particular, QD samples may be modelled 
as ensembles of two-level systems with parallel transition 
dipoles, corresponding to the fundamental optical (inter- 
band) transition between the ground ('empty dot') state 
and the confined exciton state (electron-hole pair in the 
dot), although one should not expect the transition en- 
ergies to be exactly matched in these artificial objects. 
Nonetheless, these semiconductor structures are specific 
with some respects. 

The spacing between QDs in intentionally manufac- 
tured DQD systems is typically of the order of nanome- 
ters, that is 2-3 orders of magnitude smaller than the 
wave length of the resonantly coupled radiation. This 
allows one to assume the Dicke limit of the coupling and 
neglect the retardation effects, which were one of the ma- 
jor concerns in the general theor y 25 ' 29 . On the other 
hand, small distance between the two dots precludes in- 
dividual addressing of each system by the exciting field 
so that only certain initial states may be prepared. This 
means, in particular, that the two-excitation (biexciton, 
i.e., one electron-hole pair in each QD of the DQD struc- 
ture) states must be included in the description (except 
for the weak excitation limit). 

Another important feature that distinguishes QDs 
from atomic samples is the presence of two kinds of dipole 
couplings. One of them is the direct interaction between 
static dipole moments associated with the electron-hole 
charge distribution in the two-dots^. This kind of cou- 
pling is obviously present only if the two dots are occu- 
pied by excitons and its effect is to shift the energy of the 
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FIG. 1: (Color online) Graphical representation of the energy 
levels (basis states) and couplings in the double quantum dot 
system. 



biexciton state (denoted by Vb in Fig. [T]) . The other cou- 
pling is related to the interband matrix elements of the 
electric dipole moment and is analogous to the Forster 
coupling in molecular systems as well as to dipole cou- 
plings between atoms which appear in the description of 
superradiance phenomena in atomic samples. It couples 
the two single-exciton states of the system via a process 
that may be imagined as a recombination of the exciton 
in one of the dots with the subsequent transfer of en- 
ergy to the other one (via Coulomb interaction), where 
the exciton is recreated (this coupling is denoted by V in 
Fig- 0} • Thus, this coupling has an "excitation transfer" 
character. In QDs, this interaction is modified by the 
finite size of the charge distribution (the QD size is com- 
parable with the inter-dot distance) and has a finite value 
in the formal limit of vanishing distanc o 32 i 33 ' 34 . There- 
fore, it is not a universal function of the distance. Fi- 
nally, apart from this dipole-dipole coupling, other kinds 
of interaction may be present (e.g., effective 'tunneling' 
coupling accounting for a slight overlap of wave func- 
tions which may dominate over the Forster coupling for 
closely spaced dots). For these reasons, the strength of 
the coupling between the dots becomes an essentially in- 
dependent parameter. All in all, there are three indepen- 
dent parameters governing the radiative properties of a 
double-dot system in the Dicke limit, compared to two 
in the case of atoms. 

In this paper we study the interplay of all the factors 
affecting the optical response of a DQD: the coupling be- 
tween the two dots forming the DQD, the mismatch be- 
tween their transition energies and their collective inter- 
action with the radiation environment (vacuum). With 
the recent experimental progres s 3 ^' 3 - , optical studies of 
a single, resonantly driven nanostructure have become 
feasible. Therefore, we study the simplest optical prop- 
erty of the system, that is, the response to a single ul- 
trafast pulse tuned to the fundamental interband tran- 
sition of the system. The analysis includes the linear 
regime, where the spectral properties of the optical re- 



sponse yield the absorption spectrum of the system, as 
well as the higher order contributions which are affected 
by the exciton-biexciton dynamics. We show that both 
the time-resolved signal and its frequency spectrum can 
show clear signatures of collective coupling to the radi- 
ation reservoir. Moreover, in the realistic case of non- 
identical dots, the striking features related to collective 
radiative relaxation appear only when the inter-dot cou- 
pling is strong enough and has the "excitation transfer" 
character (as opposed to the occupation-preserving biex- 
citon shift). In this way, our results provide a sensitive 
test for the appearance of a specific kind of coupling be- 
tween the dots. 

The paper is organized as follows. In Sec. HH we de- 
scribe the system under study and define its model. In 
Sec. IIII1 equations describing the evolution of the optical 
polarization are derived. Sec. IIVI contains the discussion 
of the results. Summary and final remarks are contained 
in Sec. El 



II. THE SYSTEM 

The system under study is a DQD composed of two 
QDs placed at a distance much smaller than the rele- 
vant photon wavelength. We restrict the discussion to 
the ground exciton states in the two dots. Due to strong 
electron-hole Coulomb attraction, the 'spatially direct' 
states (electron-hole pairs confined in the same dot) have 
much lower energy than the 'dissociated' states (we do 
not consider external electric fields which would change 
this pictur o 37 ' 38 ). Therefore, we include only the 'spa- 
tially direct' states in our model. We assume also that 
the polarization of the laser pulse corresponds to a po- 
larization eigenstate of the excitons, which allows us to 
include only one out of the two bright states in each dot 
and to describe the DQD as a four level system, with the 
state |0) representing empty dots, |1) and |2) denoting the 
single exciton states with an exciton localized in the first 
and in the second dot, respectively, and |3) representing 
the "molecular biexciton" state, that is, the state with 
both dots occupied by an exciton. We denote the ener- 
gies corresponding to the fundamental optical transition 
in the two dots by E12 = E ± A and allow for a cou- 
pling between the dots, whose amplitude is V. The latter 
may originate either from the interband dipole (Forster) 
coupling between the dot&i 3 - or appear as an effective de- 
scription of tunnel coupling if the carrier wave functions 
in the two dots overlap. In addition, excitons confined 
in the two dots interact via their static dipole moments, 
which shifts the energy of the biexciton state by Vb ■ The 
Hamiltonian describing the isolated DQD system is then 



iTx = (£ + A)|l>(l| + (£-A)|2><2| 

+y(|i)(2| + |2)(i|). 



(2£ + Vb)|3)<3| 



The structure of energy levels is shown in Fig. [TJ 

The two dots interact with their common radiative 
reservoir. The Hamiltonian describing this interaction 
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H int = J] d ■ e x (k) J ^~b{ x + H.c, 
with£_ = 2^ = (cri 13 +cri 2) ), where ff W = |0>(1| + |2)(3| 

(2) 

and = |0)(2| + |1)(3| are annihilation operators for 
an exciton in the first and second dot, respectively, k is 
a photon wave vector, A denotes polarizations, bk,\, b k x 
are photon annihilation and creation operators, d is the 
interband dipole moment (for simplicity equal for both 
QDs), e\ (k) is a unit polarization vector, eo is the vac- 
uum permittivity, e r is the dielectric constant of the semi- 
conductor, and v is the normalization volume for the EM 
modes. 
Finally, 



-Hphot = ^ ^kb^ x^k, 



k.X 

is the Hamiltonian of the photon reservoir, where LOk is 
the frequency of the photon with a wave vector k. 

We will describe the evolution in a "rotating basis" 
defined by the unitary transformation 

U = e iE(|lXl| + |2X2|+2|3>(3|)t/fi+ii? p hott/fi_ 

The transformed Hamiltonian is 

H = U(H X + H int + tfphot)^ + rt^U^ = H X + H inU 
where 

H X = A(|1)<1| - |2>(2|) + 7 B |3X3| + V(\l)(2\ + |2)(1|) (1) 
and 



The magnitude of the positive frequency part of the emit- 
ted coherent optical field, normalized to its initial value 
is, therefore, 

p,.s = . Pw{t) + P20(t) + P3l(t) + P32(t) (2) 
U pio(0) + p2o(0) + p 3 l(0)+p32(0)' {> 

where p is the density matrix of the four-level DQD sys- 
tem and Pij{t) = (i\p(t)\j). The matrix elements pio 
and P20 are related to the coherences between the ground 
state and the single exciton states, hence are referred to 
as exciton polarizations. The other two matrix elements, 
P3i and P32, correspond to the transition between the sin- 
gle exciton and biexciton states and are commonly called 
biexciton polarizations. 

An ultrafast pulse is assumed to be spectrally broad 
enough in order not to discriminate between the two dots 
in the structure. Due to a small (sub-wavelength) dis- 
tance between the dots, they cannot be resolved spatially, 
either. Therefore, the pulse induces optical polarizations 
symmetrically and independently in both dots. After the 
pulse, the relevant elements of the density matrix have 
the values 



Pio(O) = P2o(0) = --sin a cos 2 -, (3a) 
- — sin a sin 2 — ■ , (3b) 



P3i(0) = p 32 (0) 



where a is the pulse area. 

After this instantaneous initial excitation, the density 
matrix evolves according to the Liouville-von Neumann- 
Lindblad equation of motion 



p=-jr[H X) p]+£\p], 



with 



c[ P ] = r 



- ^{^+^-,p}- 



(4) 



(5) 



III. THE SYSTEM EVOLUTION AND THE 
OPTICAL RESPONSE 

In this section we will describe the evolution of the 
DQD system after an instantaneous excitation with an 
ultrashort pulse. An analysis of both the time-resolved 
and spectrally-resolved response will be performed. In 
the linear response limit, the latter provides the linear 
susceptibility from which the absorption spectrum can 
be derived. 

The source of the optical signal from the QD system 
is the electric dipole moment (polarization) related to 
the interband transition. Assuming identical magnitude 
and orientation of the transition dipoles in both dots, the 
relevant (interband) part of the dipole moment operator 
can be written as 

d = dY, +H.c. 



where 



r = 



E 3 \d\ : 



37reoc 3 7i 4 



is the spontaneous decay rate. It can be noted that, from 
the formal point of view, the evolution equation Q is 
not unique. In fact, a standard microscopic derivation-^ 
involves grouping the transitions into sets characterized 
by the values of the transition energies and, therefore, 
produces two different Lindblad equations, depending 
on the arbitrary classification of the two radiating sys- 
tems as "identical" or "different". Therefore, it does 
not provide any means for the study of the transition 
from the independent decay to the collective regime as 
the transition energy mismatch A is varied from zero 
to a finite value. On the other hand, Eq. ([3]), com- 
posed of the unitary and dissipative parts, is compati- 
ble with the Wigner-Weisskopf description of the double 
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QD syste m 21 ' 28 , where the Markov approximation is per- 
formed without any arbitrariness (assuming A <C E, so 
that the reservoir spectral density does not vary consid- 
erably over the relevant frequency range) . 

Note that the collective coupling of the two dots to the 
electromagnetic field, as described by Eq. ([5]), results in 
the appearance of cross-terms of the form a_ pa+ , i 7^ j, 
which are absent in the case of individual coupling to sep- 
arate reservoirs (see Appendix). These terms correspond 
to the imaginary part of the coupling that appears in 
other approache d 28 ! 30 ' 31 . Thus, in the limit of vanishing 
biexciton shift and (in the case of Refs. [30.31) of identi- 
cal dots, our model is fully equivalent to the small-system 
limit of those theories. 

The equation of motion ((4]) leads to the closed system 
of four equations for the matrix elements of interest, 



Pw = 













pw + 




-0 



P20 



+r(p 3 i + /9 32 ), 













P20 = 


K- 


~) Pro + 


(4 


-f) 



P20 



P32), 



P32 



P31 



+r(p 3 i 

Vb .A 

h 

Y. - - 

~h ~ 2 

v r 



~2 



-7T PS? 



+ 



P31, 



P32 

I -Us,. A 



3r 
T 



P31- 



(6a) 



(6b) 



(6c) 



(6d) 



sign) 

Al,2 
^3,4 
Ai, 2 

-81,2 

-43,4 

-63,4 
with 



r 

2' 



3r 

i(±n+-v B /h)- T , 



it 



v/h - iT/2 



(8a) 
(8b) 
(8c) 



T2TV B £l- + 2TV B V/h ~ iT 2 V B ~ 2A 2 T/ 



0_ [T2TV - (Wb + hT)(2n_ ± Vb/% t- «r)] ' 
0, (8e) 
±(Vb/?i - 2ir)(2iF - iV B + KT) + 2iV B n_ 



+ -~ 



(iV B + HT)(2n + T V B /h±iT) ± 2TV 

Arv 2 /h - ivgv/h 2 + rv 2 /(2h) + nr 3 



n+ [{iV B + KT) (20+ =p Vb/K ± iT) ± 2r 



n± = v/(^/?i) 2 ± ivr/& - r 2 /4 + (A/^) 2 . 

The Fourier transform of this signal is 



where 



and 



cos 2 f/(w)+sin 2 § ff (w) 
(u> — i\i)(u> — i\2)(u> — iXs)(u> — i\4) ' 



f(u>) = i(u + V/H)(lu — iX3)(u) — 1X4) 



(9) 



(10) 



(11) 



= Vb/%- V/h + 3iT)(w-i\ 1 )(<j-i\ 2 ) 

+2TV B (V/h + Lu)/h- 2A 2 T/h 2 . 

It should be noted that the frequency u> in the above 
equations (and in the following discussion) is defined with 
respect to the mean transition frequency, that is, u> — 
lu — E/h, where u> is the actual frequency of the emitted 
electromagnetic radiation. 



IV. DISCUSSION 



This system of equations with the initial values given 
by Eqs. (|3"albl can be solved by the standard Laplace 
transform method. Then, for the total emitted signal 
[Eq. ©], we find 



nt) = 5 E [ cos2 ? A * + sin2 f * 



(7) 



The evolution of the optical polarization after an in- 
stantaneous excitation depends on whether the dots are 
coupled or not. In this section we will discuss the two 
cases, comparing the optical response of the two dots in- 
teracting with the common reservoir in the Dicke limit 
(separation of the dots is small compared to Tic/E), us- 
ing the solution derived in Sec. IIII1 with the response of 
a hypothetical system consisting of two dots interacting 
with independent reservoirs. 



A. Uncoupled dots 



where (with the first index corresponding to the upper 



In the case of uncoupled QDs (V = V B = 0) the optical 
response is determined by the interplay of the other two 
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FIG. 2: (Color online) (a,c) The optical polarization as a func- 
tion of time after an ultrafast excitation for uncoupled dots 
(V = Vb = 0) in the linear response limit (a — > 0). (b,d) The 
corresponding spectrum. Figs. (a,b) show the results for QDs 
interacting with a common reservoir in the Dicke limit, while 
Figs. (c,d) refer to dots radiating into independent reservoirs. 
Here F = 0.01 ps" 1 and the values of A/h are as shown in the 
figure (in ps -1 ). Line definitions in (b,d) are the same as in 
(a,c). Dotted lines in (a,c) show the envelope ±exp(— Tt/2). 
The vertical scale in (b) and (d) is the same. 



parameters: the recombination rate T and the transition 
energy mismatch A. In Fig. [5] we show the optical po- 
larization in the time and frequency domain for a fixed 
value of the radiative recombination time 1/T = 100 ps. 
We assume here a — > (linear response limit), so that 
the imaginary part of the Fourier transform of the nor- 
malized signal is proportional to the absorption spectrum 
of the DQD. In this case, P(t) is purely imaginary. 

For A^ KT the evolution of the optical signal is dom- 
inated by optical beats due to the interference of fields 
emitted from the two dots [Fig. [2]^a,c), gray dash-dotted 
line] . There is no noticeable difference between the cases 
of a common reservoir [Fig. [2Ja)] and separate reservoirs 
[Fig. H{c)]. This is not surprising since systems with dif- 
ferent transition energies emit into different frequency 
sectors of the reservoir and thus essentially interact with 
different reservoirs anyway. The two spectra [Fig.[5fb,d)] 
also look indistinguishable in this case. 

This situation changes as the energy mismatch A is de- 
creased (Fig. [3J blue and green dashed lines). One effect, 
which follows from Eqs. (|5a|) and ([5]), is that the fre- 
quency is decreased in the collective case, although this 
does not lead to any qualitative difference in the evo- 
lution of the optical polarization represented in the time 
domain [Fig.[5[a,c)]. A much more pronounced difference 
is visible in the spectra [Fig.[^b,d)]. From Eqs. (jTUJ) and 



(TTTT) one has in the present special case, for A > hT/2 

knP(u\a = 0, V = 0, A > hT/2) = 

Tlu/Q + Tlu/Q+ 
4(w-ft+) 2 +r 2 /4 ~ A{uj + tt + ) 2 +r 2 /4' 

where 0+ = ^(A/?i) 2 - r 2 /4. This takes the form of 
two Lorentzians centered around to = ±fi+ only as long 
as r/fi_|_ is small and only for \oj — f2+| < T. In fact, 
P(lo — 0) = 0, which means that the spectrum must 
considerably differ from the sum of two Lorentzians when 
the latter overlap. This can be seen in Figs.[2Jb,d) (blue 
dashed line). 

The particular features of the DQD spectrum in the 
common reservoir case become even more pronounced 
when A < hT/2. Now fi + is imaginary and the absorp- 
tion spectrum can be written as 

lmP{uj\a = 0,V = 0,A< HT/2) = 

r + r + r_ r_ 
\n + \^ + r 2 + ~ p7| 4^ 2 + r 2 _ ' 

where r± = T ± 2|fi+|. Since T + -> 2r as A -> 0, 
the width of the first Lorentzian becomes twice larger 
than it was for independent reservoirs. This term cor- 
responds to the superradiant component of the system 
evolutional. The second Lorentzian (which is negative) 
becomes narrow, since T_ — > as A — > 0, and represents 
the subradiant component. However, its weight vanishes 
in the limit of identical dots [note that the amplitudes 
of the two Lorentzians are always exactly opposite, in 
accordance with P(uj = 0) = 0]. These spectral fea- 
tures are reflected by the evolution in the time domain 
shown in Fig. [2^a). For A slightly below KT/2 (green 
long dashed line), the evolution is a sum of two expo- 
nential factors: one positive, large, and short living and 
the other one negative, small, and long living. In fact, 
however, the resulting behavior is hardly distinguishable 
from the strongly damped oscillations appearing for in- 
dividually emitting dots [Fig. [He)]. Only for A <C KT 
the difference becomes remarkable: in the case of collec- 
tive emission, the polarization becomes dominated by the 
superradiant component which leads to a decay with a 
doubled rate, as compared to the separate reservoir case 
[red solid lines in Figs. [DJa,c)]. 

These results allow one to understand the transition 
from the limit of different systems (separate reservoirs) 
to identical systems (common reservoir). This transi- 
tion is manifested in the reconstruction of the absorption 
spectrum, as shown in Fig. [3] As long as A ^> hT/2, 
the two dots are coupled to different frequency sectors 
of the electromagnetic reservoir and the spectrum of the 
DQD system is almost indistinguishable from that of two 
dots coupled to separate reservoirs. In the most interest- 
ing parameter range, A ~ hT/2, the spectrum is non- 
Lorentzian as long as A > KT/2 and then switches to an 
unusual form of two Lorentzians with different weights 
and opposite signs, centered at zero frequency. Only 
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FIG. 3: (Color online) The transition from the 'different dots' 
regime to the 'identical dots' regime as manifested by the 
form of the absorption spectrum of a DQD system. Here V = 
0.01 ps -1 and the values of A/h are shown in the figures. Red 
solid lines: the actual system; blue dashed lines: hypothetical 
system made of two dots interacting with separate reservoirs. 
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FIG. 4: (Color online) (a) The optical polarization as a 
function of time after an ultrafast excitation for dots cou- 
pled by the static dipole moments only (V = 0, Vb ^ 0) 
for tan 2 (a/2) = 0.2, F = 0.01 ps -1 , A/h = 0.01 ps -1 , 
Vs/h — —0.3 ps -1 . Blue dashed line: real part, red solid 
line: imaginary part, (b) The Fourier transform of the po- 
larization signal. Red solid line: parameters as in (a). Green 
dashed line: A/h = 0.001 ps -1 . 

then the evolution of polarization has two exponentially 
decaying components. 

Although the features discussed above are interesting 
from a general physical point of view, their appearance 
requires that the energy mismatch between the dots is 
comparable with the radiative line width, that is, of the 
order of tens of /ieV. In spite of the rapid progress of 
nanostructure manufacturing, this can be very hard to 
achieve experimentally, unless the dots can be driven to 
resonance by external fields (e.g. taking advantage of a 
different strength of the DC Stark effect in the two dots). 
As we will show in Sec. IIV CI if the dots are coupled by 
an "excitation transfer" interaction, the collective effects 
manifest themselves already for much larger values of the 
energy mismatch. Before we proceed to this case, we 
study the effect of the other type of long-range interaction 
between confined excitons which is due to static electric 
dipoles and results in a biexciton shift. 



B. Coupled dots: Biexciton shift 

In this section we discuss the evolution of the optical 
polarization for a system of two QDs coupled by a static 



dipole interaction Vb (we keep V = 0). This coupling 
is important only for the third (and higher) order terms 
in the optical signal. In the parameter range where the 
transition between the two regimes of evolution occurs, as 
discussed in Sec. IIV Al the biexciton shift Vb (which is of 
the order of meV) is likely to be the largest energy scale of 
the problem. Therefore, we restrict the discussion to the 
case of V B > A and V B > HY. Then, from Eq. j8d|, in 
the leading order, ~ HT/Vb, which means that the 
correction to the components evolving with the frequen- 
cies iXi and 1X2, as well as to the spectrum of the optical 
response in the frequency range studied in Sec. IIV Al is 
negligible. 

In the spectral region of u> ~ Vb/K-, a new feature ap- 
pears when Vb ^ [see Fig. HJb)]. The shape of this 
feature evolves in a different way, as compared to the 
absorption line at w = 0. For A ^> TiT two nearly 
Lorentzian lines of almost equal magnitude appear at 
u) = (Vb ± A)/?i. These two lines collapse as A ~ hY but 
in this case they retain their approximately Lorentzian 
shape. For A < hY/2 both Lorentzians are centered at 
the same frequency w = Vb /h, as in the absorption spec- 
trum discussed above, but now their widths tend to F and 
2r as A — > 0. The area of the narrower line (which is 
negative) vanishes in this limit which, in view of the finite 
width, implies that the amplitude of this line vanishes. 
On the other hand, the area of the broader line remains 
finite. This means that the evolution of the biexciton 
beats [presented in Fig. 0Ja)] shows a two-rate damping 
which becomes dominated by the short-living component 
as the dots become identical. In contrast, as follows from 
the formulas listed in the Appendix, in the case of inde- 
pendent reservoir the polarization decay rate is always 
3T/2. 



C. Coupled dots: transfer interaction 

In this section we consider the case of two dots coupled 
by an "excitation transfer" interaction, that is, V ^ in 
Eq. |T]). We assume that the magnitude of this coupling 
is larger than the relaxation rate, V ^> KY. As we will 
see, values of the energy mismatch A for which the tran- 
sition to the collective behavior takes place is now of the 
order of V so that Y becomes the smallest frequency in 
the problem. One can, therefore, simplify the discussion 
by retaining only the terms up to linear order in Y in 
Eqs. (|5alb) and taking the amplitudes given by Eqs. (|8cT 
f) at T = 0. 

Then one gets A^a = — Y T /2 and A\$ = 1 =F 
V7(fifi), where T T = [l^fV/ (hn)]Y, = [A 2 + V 2 ] 1 ' 2 /h. 
Thus, in the linear response limit, the system shows a 
decay described by two exponential components. The 
subradiant one decays with the rate T_/2 < Y/2 and its 
amplitude vanishes when A <C V. The decay rate of the 
superradiant component is T + /2 > Y/2. This component 
dominates the decay in the limit of strongly coupled dots. 

The evolution of the optical polarization in the present 
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FIG. 5: (Color online) Left panels: envelope of the optical 
beats for two coupled QDs emitting to a common reservoir 
and to separate reservoirs. For each value of A/h (shown in 
the plot), the coupling V is adjusted so that Q = 3 ps -1 in 
each case. Right panels: the corresponding absorption spec- 
tra. Plots (a) and (b) show the properties of a DQD emit- 
ting to a common reservoir, while plots (c,d) correspond to 
two dots emitting to separate reservoirs. Dotted lines in (a,c) 
show exponential decay curves with the rates T/2 and F, while 
in (b,d) they show Lorentzians of equal weight and width T/2. 



case is shown in Fig.[5]Ja), where we plot the envelopes of 
the optical beats for three different sets of parameters A 
and V. Two effects that may be noticed are the decrease 
of the beating amplitude and the change in the decay 
rate of the signal as the energy mismatch decreases and 
the coupling increases. By comparison with the case of 
independent reservoirs [Fig. [5jc)] one can see that the 
first of these two effects appears in both cases. A non- 
selective excitation induces a symmetric superposition of 
the two occupations which has a larger overlap with one 
of the eigenstates of the coupled dots. Therefore, one 
of the oscillators contributing to the beats has a larger 
amplitude and the beating amplitude is reduced. 

In contrast, the other effect is related to the collective 
interaction with the photon reservoir. In the physical 
case of a common reservoir, the eigenstate dominating 
the response is the superradiant one. In the limit of 
A — > only the superradiant state is excited (it coin- 
cides with the optically active symmetric superposition) . 
This is reflected by the decay of the optical response, as 
shown in Fig.fSJa). For A ^> V, the optical polarization 
decays with the rate T/2, characteristic of a single sys- 
tem. As the coupling V becomes comparable with the 
energy mismatch A, the decay becomes non-exponential 
and, in fact, contains two components decaying with dif- 
ferent rates. When the coupling dominates over the en- 
ergy mismatch, 7> A, the subradiant component van- 
ishes and the signal decays with the doubled rate T. On 
the contrary, hypothetical dots coupled to independent 
reservoirs [Fig. |5Jc)] always show a decay of the optical 



response with the same rate T/2. 

The absorption spectrum corresponding to the 
time-resolved response discussed above is shown in 
Figs. [5]Jb,d). The two models of common and separate 
reservoirs differ essentially. In both cases there is a simi- 
lar transfer of line weight from one line to the other as the 
energy mismatch decreases and the coupling increases. 
However, the way the line shapes change is very differ- 
ent. In the case of independent reservoirs [Fig. EJd)], the 
line widths remain constant (no sub- and superradiance 
effects) and only the line amplitudes change. On the con- 
trary, when the dots are coupled to a common reservoir, 
the amplitudes of the lines are almost constant but their 
widths change. This is quite a remarkable signature of 
the joint appearance of inter-dot coupling and collective 
radiative decoherence in the system. It is also interest- 
ing to note that, since the weights of each line behaves 
almost in the same way in both cases, no difference can 
be observed if the absorption is averaged over an inho- 
mogeneous ensemble. 



V. CONCLUSIONS 

We have shown that the coexistence of coupling be- 
tween quantum dots and collective effects in their inter- 
action with the electromagnetic environment leaves clear 
traces in the optical response of these systems. 

For uncoupled dots, the collective radiative proper- 
ties (sub- and superradiance) become important only 
when the energy mismatch falls below the absorption line 
width. The decay of the polarization then evolves from 
damped oscillations (beats) for different dots to a super- 
radiant exponential decay for identical dots. For very 
similar dots the corresponding absorption line is com- 
posed of two Lorentzians superposed at the same tran- 
sition frequency with the same amplitude but different 
widths and opposite signs. 

Coupling between the dots changes this picture con- 
siderably. Now superradiance effects can be observed as 
long as the energy mismatch is smaller or comparable to 
the coupling strength. The envelope of the optical beats 
in this case decays with the usual rate for different dots 
and with the double (superradiant) rate when the cou- 
pling becomes much larger than the energy mismatch. 
In the intermediate range, the decay is composed of two 
exponential components. 

These effects in the time-resolved optical response are 
reflected in the absorption spectrum, where the lines cor- 
responding to the small and large components in the op- 
tical response loose or gain, respectively, their widths 
(reflecting the sub- and superradiance property). This 
presents a clear difference with respect to the (hypothet- 
ical) case of dots emitting to different reservoirs where the 
lines loose or gain amplitude, while their widths remain 
constant. This essentially different form of absorption 
lines which, for sufficiently strongly coupled dots, can be 
observed even for a transition energy mismatch in the 
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range of milli-electron-volts provides an experimentally 
accessible signature of coupling and collective decoher- 
ence in double quantum dots. 
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APPENDIX A: RADIATIVE DECAY TO 
INDEPENDENT RESERVOIRS 



By solving this system of equations one finds the normal- 
ized optical polarization defined by Eq. with 



r 3r 

Ai, 2 = ±ifi--, X 3A = i(±Q-V B /%)-—, 



B 



T[(V B /h - iT){Ttt + V/h) - 2A 2 /h 2 } 
2tl[(-iV B /n - T)(±V B /h TiT + 2Q)} ' 

(iv B + hr){v B - 2V)(y/n±n) + 2A 2 v B /h 

Q(iV B + hT)(±V B T i^T - 2hQ) ' 



1.2 



B3A = 



In this appendix, we derive the equations of motion for 
the optical polarizations and find the resulting optical 
response in the case of a two coupled QDs interacting 
with separate photon reservoirs. 

The initial values for the relevant matrix elements (at 
t = 0) are given by Eqs. (f3"al b). At t > 0, the density 
matrix evolves according to Eq. ([4]) but the Lindblad op- 
erator now has the form 



c\p] = £ r 



W (*) 1 r (») (*) i 



This equation corresponds to the system of evolution 
equations for the optical coherences 

Pw = (-iA/h~T/2) Pw - i(V/h) P2Q + Tp 32 , 

P20 = -i(V/h)p w + (iA/h-T/2)p 20 + Tp 31 , 

P32 = {-iV B /h-iA/h-3T/2)p 32 +i(V/h)p 31 , 

hi = i(V/h)p 32 + (-iV B /h + iA/h-3T/2)p 31 . 



where fl = \JV 2 + A 2 /h. 



The Fourier transform of this signal is 



cos 2 § f(w) + sin 2 |g(aQ 
(u> — i\i)(u> — i\ 2 )(uj — iX3)(u> — 



where f(u>) = i(u> + V/h + iT/2)(u> — i\3)(oj — 1X4) and 

g(u) = i(uj~V/h-V B /h + 5ir/2)(Lj-iXi)(u}~iX 2 ) 
+T(lu + V/h + iT/2)(V B /h - iT) - 2A 2 Y/h 2 . 
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